Model of fluorescence intermittency of single colloidal semiconductor quantum dots 

using multiple recombination centers 
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We present a new physical model resolving a long-standing mystery of the power-law distributions 
of the blinking times in single colloidal quantum dot fluorescence. The model considers the non- 
radiative relaxation of the exciton through multiple recombination centers. Each center is allowed to 
switch between two quasi-stationary states. We point out that the conventional threshold analysis 
method used to extract the exponents of the distributions for the on-times and off-times has a serious 
flaw: The qualitative properties of the distributions strongly depend on the threshold value chosen 
for separating the on and off states. Our new model explains naturally this threshold dependence, 
as well as other key experimental features of the single quantum dot fluorescence trajectories, such 
as the power- law power spectrum (1/f noise). 



Substantial progress has been made recently in the 
study of long range correlations in the fluctuations of 
the emission intensity (blinking) in single colloidal semi- 
conductor nanocrystals (QD) 0, U [3j nanorods 
0, nanowiresQ and some organic molecules By in- 
troducing an intensity threshold level to separate bright 
(on) and dark (off) states, Kuno et al. found that the 
on- and off-time distributions in QDs exhibit a spectac- 
ular power-law dependence over 5-6 orders of magnitude 
in time. 
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As discovered later by Shimizu et al. |3|, the power- 
law on-time distribution is cut off at times ranging from 
a few seconds to 100 s, depending on the dot and its 
environment. During the past eight years or so the 
truncated power-law form of the blinking on-time dis- 
tributions was confirmed by many experimental groups 
(see d, 0, H, [l^ and references therein) , but its micro- 
scopic origin remains a mystery. Remarkably, there are 
no "set" values for the on-time and off-time exponents. 
They are scattered in the region from 1.2 to 2.0. Similar 
on- and off-time distributions were found recently for the 
other blinking systems mentioned above: semiconductor 
nanorods(NRs)[a, [ll[, nanowires (NWs)[7[ and organic 
dyes 8\. The generality of the phenomenon is rather 
intriguing. We argued that there must be a common 
underlying mechanism responsible for the long time cor- 
related fluorescence intermittency detected in all these 
systems [12 1 . Most theoretical explanations of the QD 
blinking [2 . [sl. [l3 . 14, 3 are based on the Efros/Rosen 
charging mechanism [l6| . The mechanism attributes on- 
and off- states to a neutral and a charged QD, respec- 
tively. The light-induced electronic excitation in the 
charged QD is quenched by a fast Auger recombination 
process. A number of experimental results indicate, how- 
ever, that there are no unique bright (on) and dark (off) 
states of the QD, but a continuous set of emission intensi- 
ties 17 , 13, 131 • One can therefore suggest an alternative 
mechanism of the blinking, assuming slow fluctuations in 



the non-radiative recombination rate of the excited state 



In order to gain further insight into the possible blink- 
ing mechanism, we performed an extensive analysis of 
the on- and off- distributions of actual single QD fluores- 
cence trajectories. Our procedure is different from the 
conventional ones, as we applied the Maximum Likeli- 
hood Estimator (MLE) method to find the best Gamma- 
distribution p{t) ^ t^™ exp(— t/T) flt for the set of on 
and off durations. The MLE approach [2^ gives an un- 
biased estimation for the parameters of the power-law 
distribution with minimal statistical error. These prop- 
erties are crucial and allow for the investigation of a sin- 
gle trajectory. Our approach, in contrast with procedure 
used by Hoogenboom et al. [13], allows us to find opti- 
mal values for not only for m, but for the truncation time 
T as well. The fluorescence trajectories we investigated 
were obtained by Protasenko and Kuno and have already 
been analyzed by others 3, 2^. 

Our fitting procedure is performed repeatedly for a 
number of threshold values for each trajectory. In all the 
cases the off-time distribution truncation time is found 
too long to be detected. Also, the threshold depen- 
dence of the distributions was all but ignored until now. 
The only exception is the recent observation made on 
nanorods (not QDs) by Drndic and her coworkers llj . In 
any case, the fundamental nature of this dependence was 
not revealed until now. An example of the threshold de- 
pendency of the power law exponent (slope on log-scale) 
and on- truncation time for a singe QD trajectory is pre- 
sented in Fig[T] While we investigated a large number of 
trajectories, we have deliberately chosen for this paper 
one with clearly visible telegraph noise-like features and 
well-defined on and off maxima in the intensity histogram 
[see inset in Fig. 1(b)]. As it is evident from Fig[Tl even 
for this apparently ideal case, the distribution parame- 
ters are strongly threshold dependent. While the ma- 
jority of the analyzed trajectories are not like telegraph 
noise, we mention that all show similar threshold depen- 
dence. The on-time truncation time decreases monoton- 
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ically with increasing of the threshold. This trend is the 
same for most single QD fluorescence trajectories we an- 
alyzed. The scaling of the slope as a function of the 
threshold is more complicated. The exponent of the off- 
time distribution shows several extrema, whereas the on- 
time exponent has a minimum as the threshold value is 
varied. We wish to emphasize that dependence of on- and 
off-time exponents on threshold can qualitatively change 
from one trajectory to another. 
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FIG. 1; (color online) The threshold dependence of the on- 
time (red squares) and off-time (blue circles) distribution ex- 
ponents (a) and on-time distribution truncation time (b) ob- 
tained from the experimentally measured single QD fluores- 
cence trajectory. Error bars show the standard deviation val- 
ues. Insert b: a part of the trajectory and an intensity his- 
togram. 

We interpret this strong dependence on threshold as 
a clear indication that the standard trajectory analysis - 
based on the separation between on and off events with a 
somewhat arbitrary threshold - is not quite adequate, and 
the trajectories should be analyzed over the full range of 
threshold parameter. It also explains wide distribution 
of the exponents found by different groups. As shown 
below, one of the key results of this paper is to exploit 
the threshold dependence of the trajectory parameters to 
retrieve important information about the physical mech- 
anism of the fluctuations. 

The power spectrum of the fluorescence trajectory of 



a single QD has a power law form 2^, 27 1 Si{ll!) ~ w"' 
where I is close to 1. Therefore, we can consider the 
QD blinking process as an example of single particle 1/f 
(flicker) noise. The generally accepted phenomenological 
model for the electrical 1/f noise generation in solids is 
that of electrical transport in the presence of an environ- 
ment consisting of multiple stochastic two-level systems 
(TLS)[2i, liil. In the case of QD blinking we suggest a 
similar physical model based on a TLS environment [sot . 

In our model the non-radiative relaxation of the QD 
excitation occurs via trapping of holes to one of the N 
quenching centers, followed by a non-radiative recom- 
bination with the remaining electron. Each of these 
quenching centers could be dynamically switched be- 
tween inactive and active conformations. The two confor- 
mational states differ in their ability to trap holes: the 
hole trapping rate is much larger in the active confor- 
mation than it is in the inactive state. Recent studies 
of trapping rates in the single QDs [35| showed that the 
number of hole traps on the surface and on the core/shell 
interface is in order of 10. Interestingly, we find that we 
only need a similar number of recombination centers in 
order to reproduce the basic features of the fluorescence 
trajectories. A possible microscopic origin of the con- 
formation change in the recombination center could be 
due to the light-induced jumps of the surface or inter- 
face atom between two quasi-stable positions. The sur- 
face atoms in such a small object as colloidal QD can be 
found in a variety of local crystal configurations. Conse- 
quently, we can expect a wide distribution of switching 
rates. The non-radiative trapping rate in our model can 
therefore be expressed as 
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kt{t) =^fc,(T,(t) + fco 
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For each TLS the stochastic variable ai{t) randomly 
jumps between two values and 1, corresponding to inac- 
tive and active conformations, respectively. Furthermore, 
ki is the trapping rate in the active configuration, and 
ko is the background non-radiative relaxation rate. The 
time distribution functions for the = ^ 1 transitions 
and (7; = 1 — !■ transitions for the i-th TLS are expo- 
nential and can be characterized by the transition rates 
and 7~, respectively. While in the simplest model 
the transition rates for the individual TLS are constants 
(non-interacting TLSs), we will show that a more general 
case of the interacting TLS systems must also be consid- 
ered. The power spectral density of the process ^ within 
the non-interacting TLS model is a sum of Lorentzians 



1 ^ 



(3) 



The number of parameters in the above expression can 
be drastically reduced if the experimental constraint of 
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1/f noise spectrum is imposed. Indeed, after choosing 
ki = k and = j~ = = 7oa\ where a <C 1, one 
can effectively fit the spectrum in Eq. ([3]) with 1/w in 
the frequency region 7jv <C w ^ 71 [28|. Assuming low 
excitation intensities and steady-state conditions for the 
fermionic degrees of freedom, the quantum yield Y{t) is 
given bv(20j: 

where kr is the radiative relaxation rate. 

Let us now show that our suggested model of fluores- 
cence fluctuations exhibits strong threshold dependence 
of the on- and off-time distribution parameters. The 
problem of finding these distributions for the stochastic 
process Y{t) with known properties and threshold value 
y is equivalent to the well-known crossing problem [33| . 
There are only few cases when this problem can be solved 
exactly [33|. Fortunately, our present model can be re- 
duced to such an exactly solvable case. The system at any 
moment t could be completely described by the configu- 
ration E = {cti, . . . , (Tjv}. Clearly, there are 2^ different 
configurations. A random walk in the given configuration 
space is a Markovian stochastic process. The vector P 
containing probabilities of all configurations P^ satisfies 
the Master equation 



-P{t) = WP{t) 



(5) 



where the transition matrix W contains the following 
nonzero elements 



where E,^ = {cti, . . . ,ai — 0, . . . , ajy} and 'Sf — 
{(Ti, . . . , CTi — 1, . . . , (Ttv} for each given S. The non- 
radiative relaxation rate for given configuration E can be 
expressed by Eq. which allows us to find the corre- 
sponding emission intensity level Ys from Eq. ([4]). Let us 
introduce a threshold value for the quantum yield y. By 
definition, the QD is in the bright (on) state if Y{t) > y 
and dark (off) state otherwise. For each threshold level y 
all configurations can be separated to a bright group, sat- 
isfying a condition Is > y and a dark group. The vector 

of probabilities can be presented in the form ^ ~ (^p 

where vectors Pb and Pd contain probabilities of bright 
and dark configurations, respectively. The transition ma- 

Wbb Wbd 



trix can be recast in block form W = [ - - . The 

V Wdb Wdd J 

expressions for the normalized on-time and off-time dis- 
tribution functions in this notations are well-known [351 



Pon(i) = ( 1, Wdb eMWbbt)WbdPe ) ( 1, WbdPe 



where {a,b) denotes the scalar product, 1 is the unity 
vector and Pe is the equilibrium probabilities vector, sat- 
isfying a stationary condition WPe — 0. We found that 
the on-time and off-time distributions generated by Eqs. 
@ can be fitted by a power law function ([T]) (see the 
insert in Fig. 2a). Beyond a certain off/on time value 
the power law behavior sharply changes to exponential 
asymptotic behavior exp(— t/T). In our analysis, this 
value of T is defined as a truncation time. We performed 
simulations of on-time and off-time distributions for the 
model system of non-interacting TLS. This relatively 
simple model reproduces the general trend seen exper- 
imentally in the truncation times: The on-time trunca- 
tion decreases and the off-time truncation increases when 
the threshold value goes up. 

While this simple, non-interacting TLS model is use- 
ful in illustrating our procedure, it cannot reproduce the 
threshold dependence of the exponents. The slope of 
the on-time distribution monotonically increases with the 
threshold value, when the off-time exponent has an op- 
posite trend. In order to make our model more realistic, 
we introduce interaction between TLS in the simplest 
mean- field form (similar to Ref.[3l|). The interaction is 
characterized by the parameter a, whereas the bias for 
an individual TLS is parameterized by /3: 
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Fig. 2 provides the numerical calculation results for 
the interacting model with the following parameters: 
iV 10, 71 = 1, a = 10-1/2, kr/k = 1, fco = 0, a = 0.27 
and P = —0.13. As seen from this figure, the thresh- 
old dependence of the truncation times keep the same 
trend as for noninteracting case. In contrast, the slopes 
now show a non-monotonic threshold dependence repro- 
ducing qualitatively the experimental behavior shown in 
Fig. la. The insert in Fig. 2b shows that the interacting 
TLS model is capable of generating the two-maximum 
intensity distribution seen in Fig lb. The relative ease 
with which our simple phenomenological model captured 
the experimental trend gives us hope that the model can 
be used to extract interaction parameters for the TLS 
environment. These parameters could provide useful ex- 
perimental constraints on future microscopic models for 
the TLS environment of a variety of systems showing 
fluorescence intermittency. The model proposed here 
also explains recent observations of the non-blinking dots. 
Furthermore, a similar model can be constructed for the 
fluorescence intermittency seen in quantum wires. The 
details for these results will be published elsewhere. 

In conclusion, the phenomenological model we pro- 
posed in this paper succeeds in qualitatively explaining 
the key experimental facts characterizing long-correlated 
fluorescence intensity fluctuations of the single colloidal 
quantum dots: (1) the truncated power-law distributions 
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FIG. 2; (color online) The theoretical threshold dependence of 
the on-time (red squares) and ofF-time (blue circles) distribu- 
tion exponents (a) and truncation times (b) for the interact- 
ing TLS model (0 . Insert a: the on- and off-time distribution 
functions at the threshold value y=0.25. Insert b: probability 
distribution function (PDF) of the quantum yield. 



for on- and ofF-times obtained by the commonly used 
threshold procedure; (2) the strong threshold dependence 
of the distribution parameters m and T and wide range 
of the the extracted exponents; (3) the 1/f noise form of 
the power spectrum of the intensity fluctuations; (4) the 
continuous distribution of emission intensities and exci- 
tation lifetimes; (5) the weak temperature dependence 
of the fluorescence intermittency due to the light-driven 
character of the TLS switching process. 
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